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A (1+1) -dimensional quantum field theory with a degenerate vacuum (in infinite vol- 
ume) can contain particles, known as kinks, which interpolate between different vacua and 
>• . have nontrivial restrictions on their multi-particle Hilbert space. Assuming such a theory 
to be integrable, we show how to calculate the multi-kink energy levels in finite volume 
given its factorizable S'-matrix. In massive theories this can be done exactly up to contri- 
butions due to off-shell and tunneling effects that fall off exponentially with volume. As 
a first application we compare our analytical predictions for the kink scattering theories 
conjectured to describe the subleading thermal and magnetic perturbations of the tricrit- 
ical Ising model with numerical results from the truncated conformal space approach. In 
particular, for the subleading magnetic perturbation our results allow us to decide between 
the two different S'-matrices proposed by Smirnov and Zamolodchikov. 
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1. Introduction 

A unique property of 2-dimensional Minkowski space is that a double light cone divides 
the set of all space-like points into two disconnected components in a Lorentz-invariant 
way. In a (l+l)-dimensional quantum field theory (QFT) with degenerate vacuum it is 
therefore possible that the expectation value of a field approaches different values at the 
two space-like infinities. If the corresponding excitation creates a stable particle in the 
theory it is generically known as a "kink"; it is labelled by the two (possibly identical) 
vacua between which it interpolates, in addition to its energy-momentum and perhaps 
other quantum numbers. 

Some novel and interesting features arise when one considers the multi-particle Hilbert 
space of a theory with kinks. Since a multi-kink state corresponds to a sequence of vacua 
on a line, one cannot arbitrarily compose single kinks to obtain an allowed state. In other 
words, there are restrictions on the multi-kink Hilbert space — it is not a (direct sum of) 
bosonic or fermionic Fock spaces. 

Note that the solitons of the sine-Gordon model are not quite examples of what we 
want to call a kink. In this theory field configurations differing by a period of the cosine 
potential are identified, so that there is only one type of soliton and anti-soliton with no 
restrictions (except the exclusion principle) on their multi-particle states. (See however 
sect. 5 for a discussion of some variants of the sine-Gordon model where such restrictions do 
arise.) We reserve the term kink for excitations with nontrivial restrictions on their multi- 
particle states. Recently several @ @ 00 factorizable scattering theories of massive 
kinks have been proposed to describe certain integrable perturbations of conformal field 
theories (CFTs). 

The purpose of this paper is to analyze the multi-kink spectrum of such theories in 
finite volume with periodic boundary conditions, i.e. on a circle of circumference R, say. 
We will show how to calculate the energy levels exactly — up to terms due to tunneling 
and off-shell effects that are exponentially small for large R — in terms of the factorizable 
S'-matrix of the theory. 

A theoretical understanding of the finite-volume spectrum of a QFT provides a handle 
on the dynamics of a theory, allowing one, for instance, to test a conjectured exact S- 
matrix. The idea is to compare the analytical large-volume predictions obtained from 
the S'-matrix with one of the available numerical methods of calculating the small- volume 
spectrum of a QFT, e.g. lattice simulations or the "truncated conformal space approach" 
(TCSA) for a perturbed CFT (cf. subsect. 3.3). 
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There are other methods of checking a conjectured ^-matrix, involving comparison 
with the finite-volume energy levels of 0- and 1-particle states. For a generic theory of 
kinks these tests are not so easy to perform, and here the multi-particle spectrum can 
provide a very simple way of checking a proposed S'-matrix. In most of the remainder 
of this introduction we will briefly review the basic physical phenomena underlying the 
finite- volume spectrum of a QFT, and what information about the S'-matrix of the theory 
can be gained by analyzing different parts of the spectrum. 

There are basically three physical effects contributing to the finite- volume energy levels 
of a QFT, which we will refer to as (i) tunneling, (ii) off-shell, and (iii) boundary and 
scattering effects. 

(i) : Tunneling, well-known from quantum mechanics, occurs in theories with a degenerate 
vacuum in infinite volume. Since the 'energy barrier' between the vacua is proportional 
to the volume of the world, there will be instantons in finite volume, tunneling between 
the vacua and lifting their exact degeneracy. In simple situations, like </> 4 theory in d + 1 
dimensions with 'space' a hypercube of volume R d , the splitting (in the phase of spon- 
taneously broken symmetry) is of 0(R~ 1 R d / 2 e~ aR ). Here the exponential is simply 
the Boltzmann factor of the instanton action, the power R d l 2 is due to zero modes, and 
-R -1 is due to 1-loop fluctuations. In 1+1 dimensions the exponential is e~ mR , m being 
the mass of the kink. The prefactor in the above formula can be calculated perturbatively 
(see the last ref. in |J) in a given QFT, but it is apparently not known how to express 
it in a "universal" manner in terms of the S'-matrix of a generic theory — in fact, it is 
not clear if such a universal formula exists. If there is tunneling, it is not only relevant 
for the energy levels of 0-particle states; its effects on other levels will be of comparable 
magnitude, although we are not aware of any explicit calculations. 

(ii) : There are two ways in which off-shell effects influence the finite- volume spectrum. 
The first, relevant for multi-particle states, reflects the fact that in finite volume inter- 
actions between physical particles (due to the exchange of virtual particles) cannot be 
completely expressed in terms of .S-matrix elements, which connect asymptotic states. 
Since in a massive theory interactions decrease exponentially with the distance between 
the particles, the contributions of this effect to the energy levels of stationary scattering 
states (cf. (iii) and subsect. 2.2) are expected to decrease exponentially with volume. The 
second off-shell effect is due to the fact that in finite volume virtual particles can "travel 
once or several times around the world" before annihilating again or being absorbed by a 
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real particle. The resulting volume dependence of these vacuum polarization effects again 
gives rise to exponentially small corrections to the energy levels. For the case of zero mo- 
mentum 1-particle energies, i.e. finite- volume masses, these corrections have been studied 
in some detail |10[ [11]] for essentially arbitrary massive QFTs (not only integrable theories 
in 1+1 dimensions) without a degenerate vacuum in infinite volume. The leading terms 
contributing to finite-size mass corrections in such theories are given by a universal formula 
involving only S'-matrix elements. It is not known if the exact finite-size mass corrections 
can also be expressed solely through scattering amplitudes. 

For O-particle states stronger results can be obtained in certain cases. In 1+1 dimen- 
sions there is a beautiful method, known as the thermodynamic Bethe Ansatz [|12| [[L3| [[Uj |T5j] | 



(TBA), to calculate the exact finite- volume ground state energy Eq(R) of an integrable 
QFT in terms of its ^-matrix. (And with some modifications this method can also be ap- 



plied to certain O-particle states above the ground state |^6[||17||HI8[| .) Although restricted 
to integrable theories, this method does not require the existence of a unique vacuum. The 
exact calculation of the ground state energy on a circle is possible due to its relation to the 
free energy on an infinite line at finite temperature (and zero chemical potentials). The 
final result for Eo(R), expressed in terms of the solution of a set of coupled non-linear 
integral equations, can provide a very strong test of a conjectured S'-matrix, since the 
small- R behaviour of Eq(R) allows one to extract the central charge and other properties 
of the CFT describing the UV limit of the QFT under consideration.0 Unfortunately, if 
the S'-matrix of the theory is not diagonal, the TBA requires a detailed case by case anal- 



ysis (see [[T9| []2(J for the first examples) . Such an analysis can be quite nontrivial for kink 



theories, and therefore other methods of checking a conjectured S'-matrix are desirable. 

(iii): From quantum mechanics one is familiar with the fact that boundary conditions 
lead to the quantization of energies and momenta. For 1-particle states in a QFT this 
effect, together with contributions due to tunneling and vacuum polarization, completely 
determines the energies in finite volume. For multi-particle states one must also take 
scattering effects into account, which simply refer to the fact that energy eigenstates in 



1 However, in some cases Eq(R) does not uniquely determine the UV CFT; this happens 
e.g. for certain perturbations of minimal CFTs with the same central charge but different modular 
invariant partition function. 
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finite volume are stationary scattering states.! The resulting quantization conditions, 
which we will refer to as Bethe Ansatz equations, can be expressed solely in terms of the 
S-matrix of the (factorizable) theory if we neglect the nonzero range of the interactions. 
The corresponding finite-size corrections to energies are generically of 0(1/ R 2 ), as we will 
see, and therefore much larger than the exponentially small (in massive theories, at least) 
contributions due to the finite interaction range, off-shell, or tunneling effects. Multi- 
particle energy levels are therefore ideally suited to provide a check on kink S- matrices. 

Previously multi-particle states were only studied for theories with diagonal S- 
matrices ||q| ||22|l ; 2-particle levels were also discussed for the 0(3) non-linear cr-model [23 



where the S'-matrix is diagonal in a 2-particle basis of isospin eigenstates. 

The paper is organized as follows. In sect. 2 we describe how to obtain the energy 
levels and eigenfunctions of multi-kink states in finite volume (up to the exponentially 
small corrections discussed above). Sects. 3 and 4 are devoted to examples. In sect. 3 
we start with a description of some general features of ^3-perturbed unitary minimal 
CFTs, leading (in one direction) to integrable theories of massive kinks. As a warm-up, we 
discuss in subsect. 3.1 the simplest member of this family, the Ising field theory in the phase 
of spontaneously broken symmetry. For this theory the complete finite-volume spectrum 
is known exactly, and provides an instructive check on the large-volume Bethe Ansatz 
predictions (which are rather trivial to obtain in this case). Subsect. 3.2 treats the next 
theory in the family, the subleading thermal perturbation of the tricritical Ising CFT. In 
subsect. 3.3 we compare our analytical predictions for 2- and 4-kink states with numerical 
results of the truncated conformal space approach (TCSA). Sect. 4 deals with another 
family of integrable theories of kinks, the 02,i-perturbed unitary minimal CFTs. We focus 
our attention on the subleading magnetic perturbation of the tricritical Ising CFT, for 
which two different S'-matrix theories have been proposed, one by Smirnov || and the 
other by Zamolodchikov . We clarify the relation between the two S-matrices and point 
out some a priori problems with Zamolodchikov's conjecture. Comparing our large-volume 
results for 2- and 3-kink energy levels with TCSA data provides very strong support for 



2 In integrable theories states with a definite number of particles will not mix with states of 
different particle number, whereas for non-integrable theories this will be true generically (i.e. in 
the absence of appropriate selection rules) only for 2-particle states with energies below the 3- 
particle threshold. See for a detailed study of 2-particle levels in higher-dimensional (and 
therefore non-integrable) theories. 
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Smirnov's proposal (although the subtle issue of "crossing factors" , cf. subsects. 2.3 and 4.1, 
remains open). In sect. 5 we present our conclusions, and discuss some simple modifications 
of the sine-Gordon model which lead to theories with kinks. 

2. The iV-Particle Transfer Matrix and the Bethe-Yang Equations 

2.1. S -matrix conventions 

In an integrable massive (1 + 1) -dimensional QFT all scattering amplitudes factorize 
into 2-particle amplitudes S^f(9 = |6>i 2 |) describing scattering processes a{9\) + b(9 2 ) — ► 
0(6*2) + d(9i) among the different particles a, 6, ... of the theory (see fig. 1(a)). 

* x 



Fig. 1: (a) The 2-particle scattering process a{6\) + b{9 2 ) — > c(0 2 ) + d{6\) 
(b) The 2-kink process K ai (9 1 ) + K lf3 (9 2 ) -> K a5 {9 2 ) + K 5j3 {9 1 ) 

Here the Oi are rapidities, parametrizing the energy-momentum of a particle on its 
mass shell via 

ip iP 1 ) = (wcosh^, msinh^) , (2.1) 

and 9ij = 9i — 9j. 

We will be particularly interested in scattering theories of kinks. For such theories it 
is more convenient to label the S'-matrix elements by the asymptotic vacua between which 
the kinks interpolate. (To be sure, by essentially just rewriting our results given below 
in kink language into "ordinary" S'-matrix notation, they apply to arbitrary factorizable 
scattering theories.) We will consider the case of a single multiplet of kinks of mass m > 
which are all bosons or all fermions, in order to keep the notation simple and transparent; 
this assumption can be trivially relaxed. The kink interpolating between the vacuum a at 
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x — > — oo and (3 at x — > +00 will be denoted by K a g(6)M We let S^piOw) describe the 
process K ai {6{) + Kypfa) — » -^05(^2) + Ksp(0i), where 6>i > #2- The latter restriction 
ensures the consistent ordering a, 7, j3 (a, S, 0) of the vacua along space in the far past 
(future), as is clear from the graphical description of the scattering process presented in 
fig. 1(b). (It will be convenient below to let a, (3, 7, 5 range over all vacua and set S^i(0) = 
if one or more of the kinks in fig. 1(b) does not exist in the theory.) 

To get used to the kink notation it is instructive to rewrite the standard require- 
ments of S'-matrix theory in this language. We will assume our theories to be time- 
reversal and parity invariant, i.e. S^giO) = S^ g (6) = Sl^(9). Using real analyticity 
^o$(0) = i^cfpi^ 9*))*, unitarity can then be written as 

E3I?(0)^(-0) = V for all Q) A 7) i. (2.2) 

The requirement of crossing symmetry amounts to 

Sl & p{0) = S${iv-6) . (2.3) 

The factorization or Yang-Baxter equations [|4j , expressing the consistent factorization of 
an arbitrary amplitude into 2-particle amplitudes, read 

5>3?(*ia) S2(0i8) S^(6 23 ) = 5>^(0 23 ) S^(d 13 ) S%(d 12 ) , (2.4) 
which is illustrated in fig. 2. 




Fig. 2: The Yang-Baxter equation (2.5) (the sum over \i on both sides of the equation is 
suppressed) . 

3 We will always use greek indices for the asymptotic vacua; this allows one to distinguish 
"ordinary" S-matrix elements, which will be adorned by latin indices for the particles, from those 
in kink notation. 
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Finally, if the kinks in a theory form (stable) bound states, as indicated by simple 
poles (with residues of appropriate |^S[ sign) of the scattering amplitudes in the bound 
state region 6 G (0,z7r), then the S'-matrix also has to satisfy bootstrap constraints pjfl fl7|. 
In theories of a single multiplet of kinks this can of course happen only if the bound states 
of the "constituent" kinks are again kinks in the same multiplet (see sect. 4 for an example). 
We will treat here only this case; generalizations are straightforward. 




Fig. 3: The bootstrap equation (2.5) is obtained by taking the residues of both sides of 
the depicted "equality" at Q\ — 9 2 = 2jri/3 (the sum over \i on the lhs is suppressed). 



Let K ai be a bound state of K a p and K Pl (in the direct channel). Considering the 
3-kink process K a p(0i) + K Pl {d 2 ) + K lS (0 3 ) -> K arj (d 3 ) + K m ($ 2 ) + K eS (9i) with 6 X and 
02 tuned to the complex values necessary to create the intermediate bound-state kink K a p 
on shell (see fig. 3), one concludes that the following equations must hold: 



£ 9^ 9am S%[0 + j ) S%{6 - j ) = g afh g V£S S%(0) (2.5) 

(where g a f3j can be cancelled on both sides whenever it is nonzero). Here the couplings 
g a /3y are defined (at worst up to signs) in terms of the residues at the bound-state poles, 

Res S^(9) = i g aPl g a ^ . (2.6) 

#=27ri/3 

Note that (|2.3| ) implies that the residue of the crossed-channel pole is just the opposite of 
the direct- channel one ( p.6| ), 

Res SV(6) = -i g aPl g aflJ . (2.7) 

8=ni/3 

If all couplings (7 Q/ 3 7 are real, the S'-matrix is said to be "1-particle unitary". Parity 
invariance, crossing symmetry, and the bootstrap equations imply that the (g a p-y) 2 are 



7 



symmetric in a, f3 and 7, as long as all S^(0) = — 1 (cf. for an analogous statement 
in diagonal S'-matrix theories of ordinary particles). The latter presumably holds in all 
bosonic kink theories, see eq. (2.10) below. For fermionic theories, where S^(0) = +1, 
there are some subtle signs that have to be taken into account, already in the definition of 
the couplings ( |2.6|) (cf. |25j for details in the case of "ordinary" scattering theories). 

2.2. The Bethe-Yang Equations 

Consider a factorizable scattering theory on a finite space with periodic boundary 
conditions, in other words on a cylinder. We are interested in multi-particle energy levels 
as a function of the "volume of the world", i.e. the circumference R of the cylinder. In 
an integrable theory particle number is conserved and so we can talk about iV-particle 
states for any fixed iV — at least as long as the notion of particles makes any sense in 
finite volume, namely when R ^> 1/m and the particles are far apart. One should think 
of an iV-particle energy eigenstate in large volume as a stationary scattering state, i.e. a 
superposition of iV-kink states in our case, which is invariant under the scattering of the 
various kinks. Such a state is characterized by a set of rapidities {6q, . . . , 0/v}, since they 
are conserved in the scattering processes. (Note that because of "mixing" between different 
particles of the same mass, the 9i cannot be assigned to specific particles). 

Let tjj = ijj(xi, . . . , xn) be the wavefunction of an iV-particle state, and let ifj a = 
% j ) a. 1 ...a, N denote its components with respect to the basis states \K aiOL2 (9i)K a2a3 (62) ■■ ■ K aNOtl (6>at)).| 
Thinking of the basis vectors as zn-states we must have 6\ > 62 > ■ ■ ■ > On- The number 
of iV-particle states (for fixed Qi) with periodic boundary conditions will be denoted by 
oIn- 

The behaviour of ifi when its arguments Xi are far apart can be expressed in terms 
of scattering amplitudes. Periodic boundary conditions are imposed by, roughly speaking, 
stipulating that the wavefunction change only by a factor of (— 1) F when "carrying a kink 
K 0Lk0lk+1 {6k) once around the cylinder" .0 (Here (— 1) F = ±1 for kinks that are bosons or 
fermions, respectively.) This leads to the following equations for the allowed if) 

e iRmsinhe * T k {6 1 ,...,0 N )Pl> ' = (~1) F ^, k = l,...,N, (2.8) 

a: a k =p k+1 

Of course, this should not be literally understood as a physical process, but rather as a 
mnemonic describing a sequence of rewritings of the wavefunction. 
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where 

N 

f k {e x ,...,e N t = n s££S (**-*<) • (2-9) 

Here a/v+i = 0,1 and /3at+i = Pi- (A graphical illustration of a matrix very closely related 
to Tk will be given below.) Note that there is no factor corresponding to i = k in ( |2.9| ), 
since u K akak+1 does not scatter with itself when taken around the circle". We will refer 



to (|2.8|) as the Bethe- Yang equations (cf. |27j] for analogous equations in a non-relativistic 
theory of "ordinary" particles). 



We would like to rewrite ( |2.8| ) in a more convenient way, getting rid of the restriction on 
the sum over the multi-index a. To do so, note that unitarity implies Ylsi^Z^i^)) 2 = 1- ^ n 
all kink theories we know this is satisfied in the simplest possible way, namely the S'-matrix 
at zero relative rapidity satisfies 5^(0) = 5^(0) 5 jS = ±5 lS . [In ordinary S'-matrix 
notation this reads £^(0) = ±<5 ac <5 M , total reflection in the low-energy limit.] 

Furthermore, as a slight generalization of the arguments and observations in [|nj , 



we claim that all (non- vanishing) amplitudes for the scattering of kinks in the same mul- 
tiplet satisfy 

S^O) = ~(-l) F , (2.10) 

since it amounts to an exclusion principle, i.e. 6i 7^ 9j for all i 7^ j. Eq. ( |2.10|) and the 
analogous statement in non-kink language is true in all consistent (l + l)-dimensional QFTs 
we are aware of, except for free bosons, which are somewhat singular from various points 
of view. 

We can now rewrite ( ^78|) as 

e iRmd*hO h T{d k \d U ...,B N f a ^ a = fc=l,...,JV, (2.11) 

a 

where there is no restriction on the sum over a. Here T is the N-particle transfer matrix, 
the analog of the (inhomogeneous) row-to-row transfer matrix (of an IRF model in our 



kink context) in statistical mechanics |28 



TV 



T{et = T(e\9i,...,e N )P = J] s£*£(0-0 4 ) , (2.12) 

presented graphically in fig. 4. 
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Fig. 4: Graphical description of the iV-kink transfer matrix (2.10). Each vertex corre- 
sponds to a 2-kink scattering amplitude, with the time axis to be thought of as pointing 
in the NE-direction. 



When calculating the finite- volume ground state energy of the theory using the ther- 
modynamic Bethe Ansatz technique, it is the iV — > oo limit of this transfer matrix that 
must be considered. Taking the thermodynamic limit, where also R — > oo with N/R 
fixed, certain simplifications occur which often allow one to determine the form of the 
dominant (in the thermodynamic sense) eigenvalues, and derive integral equations for the 



distributions of the 6i minimizing the free energy; see Jl9] [2(J for details and examples. 



Returning to finite N, for given R the Bethe- Yang equations (|2.11| ) have solutions 
only for special 9k = 9k(R). These equations therefore determine the wavefunctions ip and 
the (total) energies and momenta 



N N 



E = mcosh#fc , P = msinh^fc (2-13) 

k=l k=l 

of finite-volume eigenstates. 

A comment on the range of validity of eq. ( |2.13| ) is necessary. The Yang-Baxter 
equations and unitarity imply that for any permutation Q 

f Ql (e)-f Q2 (6).....f QN (6) = 1 , (2.14) 
so that eqs. (|2.8f) lead to the quantization of the momentum P of a state as given in ( |2.13|) , 

e iPR = ((_!)**)* . (2.15) 

This equation is in fact exact in any QFT, irrespective of the Bethe- Yang equations, 
because of translational invariance and periodicity. Therefore eq. ( |2.13[ ) for the momentum 
is exact. On the other hand, the above expression for the energy of a state is of course not 
exact for finite R — we are neglecting contributions due to the fact that the behaviour 
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of the wave function ip can be expressed in terms of scattering amplitudes only when all 
particles are far apart, as well as tunneling and polarization effects (cf. sect. 1). Fortunately, 
these effects decay exponentially with R in a massive theory, whereas for a multi-particle 
state the leading finite-size corrections to the energy levels determined by (|2.11|) - (|2.13| ) are 
0(1/ R 2 ), see subsect. 2.4. All powerlike dependence on R is contained in ( p.ll|) - (|2.13|) . 

2.3. General properties of T (6) 



The most important property of the transfer matrices T(0), eq. ( |2.12|) , is that they 



commute for different 0, with fixed 0±, . . . ,0n- This follows from the Yang-Baxter equa- 
tions. Any commuting family of transfer matrices, such as that occurring in ( |2.11| ), there- 
fore has a common set of eigenvectors tp^ s \0i, . . . , 9n) with eigenvalues \( s '(0\0i, . . . , On), 
s labelling the different eigenvectors.il In the cases we will consider the transfer matrix 
will be meromorphic in all its arguments (with no poles for real 0, Qi). 

Real analyticity and unitarity of S(0) imply that the transfer matrix T(0) is unitary 
if all — Oi G R. The solutions of (|2.11| ) we are interested in correspond to physical states, 
so they involve only real Oi. Hence all eigenvalues of T in ( |2.11| ) are of magnitude 1. This 
is of course necessary for self-consistency of Q2.11|) . 



Let us denote {6q, . . . , On} by 0, and introduce the notation 

T k (0) = T(0 k \0 l ,...,0 N ) 
\[ s) (0) ee X^(0 k \0 1 ,...,0 N ) (2-16) 
4 S) (0) = 6M(6 k \6 1 ,...,6 N ) = -ilnXi s \0) 



5 In general it is of course not easy to diagonalize the transfer matrices (note that the S-matrices 
and therefore the transfer matrices of kink theories are never diagonal). In theories with a large 
global symmetry, in particular non-kink theories, this task simplifies. Consider, for instance, the 
0(3) non-linear cr-model discussed in [23], where the spectrum consists of an isospin 1=1 triplet. 



In the 9-dimensional space of 2-particle states we can change to a basis of eigenstates of / and 
^3, in which the 2-particle transfer matrix is diagonal. Furthermore, since the basis change is 
accomplished by a matrix with rapidity-independent coefficients, the eigenvalues (which are I-, 
3- and 5-fold degenerate) will simply be linear combinations of the 2-particle S-matrix elements. 
In contrast, for a generic factorizable S-matrix theory there are not enough global charges to 
diagonalize even the 2-particle transfer matrix, and consequently the change to a diagonal basis 
will involve (9-dependent coefficients. 
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for the transfer matrix, its eigenvalues and their phases, respectively. The 27rZ ambiguity 
in the phases has to be fixed at some, say real, 9 by some convention; everywhere else the 
phases are then uniquely determined. 

Next note that real analyticty implies 

T k {-9) = T k (9)* for 9eK N . (2.17) 

Therefore the eigenvalues come in pairs X k (9), X k (9) related by 9 —>■ —9 and complex 
conjugation, or they are invariant under these operations (we then set s = s) 

\ { k s \-9) = \<p{0)*, si s) (-9) = -SM{0) (mod27r) for 9 G K N . (2.18) 

Here s = s is possible only if A^^O) G {±1}, obviously. That the converse is also true 
follows from Schwarz's reflection principle]! The same argument also shows that if we 
define 

~di s \9) = di s \9)-d^(0) , (2.19) 
where (^^(O) = 5 k s \o), which is clearly independent of /c, then 

6£>{0) = ~5 { k s \9) = -~5i S \-0) for 9eK N . (2.20) 

Let ^(9) be a solution of ( ^TT|) , T k {9)^{9) = \ { k s) (9)i(;(9), k = 1,...,N. Complex 
conjugation, using (|2.17| )-( p.l8|) , then gives for real 9 

T k (-$)i,(ff)* = \P{-0M0)* , k = l,...,N. (2.21) 

We thus see explicitly that solutions to the Bethe-Yang equations come in pairs related by 
parity, as expected. Namely, if 9 is a solution so is —9. These two solutions are distinct if 9 
and —9 are different as unordered sets; they are exactly degenerate in energy (this will also 
be true for the exact levels when tunneling and off-shell effects are taken into account) and 
belong to sectors of opposite total momentum. This is useful, for instance, in identifying 
degenerate and non-degenerate levels in the zero momentum sector. To give an example 
relevant for sect. 4.2, note that in this sector 3-particle states with rapidities {9, 0, — 9} will 

6 To prove this we need to know that the eigenvalues X < jf\6) are real for at least part of 
the imaginary 9 "axis". In the cases of physical interest, where solutions to the Bethe-Yang 
equations ( |2.8| ) exist for real 9, and by analyticity also in a complex neighborhood of 9 = 0, it is 
clear from ( p.8| ) that the X k (6) have this property. 
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not be degenerate, whereas a level {6>i, #2, 6*3} with all 9% 7^ will be degenerate with the 
state with rapidities {—6*3, —6*2, — 6*i}.0 

There is one more general property of the transfer matrix that we should mention. 
As we will see in sects. 3-4, it is of some physical interest to consider what happens if we 
change the S'-matrix as follows: 

Here the Xa{9) are functions which are, say, meromorphic in 9, but otherwise arbitrary at 
this point. From a lattice model point of view, where the S'-matrix elements correspond 
to Boltzmann weights of an IRF model, such a change might be referred to as a "gauge 
transformation" , since the partition function with periodic boundary conditions does not 
change, obviously. A slightly stronger statement is that the eigenvalues of the transfer 
matrix T{9) do not change under such a transformation (the eigenvectors do, though). 
This is physically interesting since it implies that the finite-volume multi-kink spectrum 



of two scattering theories differing just by factors as in ( |2.22| ) will be identical, though 



perhaps only up to terms that are exponentially small for large volumeJl 

To explicitly prove that the eigenvalues of T(9) do not change, note that under the 
above transformation 

_^ TT Xa.i(9 — 9i) XP i + i(& — &i) rp (a \B (0 O o\ 

fj[ Xfc (9 ~ 9i) Xa i+1 [9 - 9i) 

This implies that the diagonal elements of any power of T(9) do not change. In particular, 
the trace of any power of T(9) stays the same. Since the traces of the first n powers of an 



7 There is another mechanism by which exactly degnerate levels can arise, namely if the com- 
mon eigenspace of all Tk(6) is 2- or higher- dimensional for some eigenvalue Aj|, (0) (at least on 
some sufficiently large submanifold of 6 £ R w ). This happens in theories with a sufficiently large 
global symmetry, cf. footnote 5. 

8 And using the results of [|i~0f [p"T|| one can check that the same is true for the leading large- 
volume finite-size mass corrections (ignoring tunneling effects, for which we are unable to make 
such a statement, cf. sect. 1). 
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n x n matrix provide a complete set of invariants Jl we have proved that the eigenvalues of 
T{9) do not change. 

In physical applications the interest in ( |2.22|) is due to the following. Methods like the 
quantum-group approach (cf. sect. 3 for some brief remarks) allow one to conjecture "proto- 
S'-matrices" for certain QFTs. These objects are real analytic and satisfy the Yang-Baxter 
equations, but are not necessarily unitary or crossing symmetric. Unitarity can often be 
restored by multiplying the "proto-S'-matrix" by a suitably chosen function of 9. The 
crossing properties can be changed by factors as in ( |2.22| ). For the transformed S- matrix 
to still satisfy the Yang-Baxter equations and be real analytic, the Xa{9) have to be of the 
form Xa(0) = p a ^ 2lT ^ with p a > (up to an a-independent function of 9 which cancels 



in ( |2.22p and can therefore be ignored) . For \a ($) °f this form the crossing properties are 
affected as follows 

s .. //>„/>.; '"V../ 1 / 2 24 s 



S$(i*-0) V PyPs S$(i*-6) ' 

Note that for theories in which the kinks are bound states of themselves the couplings 



defined in (|2.6|) transform as g a ^ — > ( Pa z~' Y^gap-y under ( |2.22|) , assuming the above form 



of the "crossing factors" Xa(9). In particular, the symmetry of the couplings in a, {3, 7, 
and the equality (up to a sign) of the direct and crossed channel residues of bound-state 
poles, are valid only if strict crossing symmetry ( |2.3| ) holds. On the other hand, one can 
check that the bootstrap equations (|2.5| ) are not affected, given the above Xa(9). Examples 
where the crossing factors are important will be encountered in sects. 3-4. 

2.4- Solution of the Bethe Ansatz equations 

For any given 6 = {9±, . . . ,0n} the iV transfer matrices Tk(0) will have djsr common 
eigenvectors ^ s \6) with eigenvalues {0) of magnitude 1. For a given R and eigenvector 
ip( s \6) the Bethe- Yang equations restrict the allowed 9k of a physical state to the solutions 
of 

rsinh# fc + S^(e k \e 1} ...,e N ) = 2nn k , n k G Z + \ , k=l,...,N, (2.25) 



9 Except possibly for a discrete set of exceptional values of the traces. But in our case, where 
the eigenvalues are (branches of) meromorphic functions, the possibility of such exceptional values 
can be ignored. 
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which we will refer to as the Bethe Ansatz (BA) equations, to distinguish them from the 
Bethe-Yang equations ( |2.11| ) . Here we introduced the dimensionless "volume" of the world 
r = mR, measured in units of the Compton wavelength of the kinks. 

We assume that for given s and {n k } eqs. (|2.25[ ) have a unique solution 9% = 9 k (r), k = 
1, . . . , N. At least physically this is clear, by "continuity of the infinite-volume limit": As 
r — > oo all 9 k — > 0, and the unique solution is 9 k = (2nn k - 5 (s) (0))/r + C(l/r 2 ). 

The level of "type" s and "quantum numbers" n k will be denoted by s(ni, . . . ,njv). 
The allowed ni, . . . , n k E Z + | are subject to the constraints that in a sector of total 
momentum P = 2nn/R 



k=i 

and ni > ri2 > . . . > njvj to comply with the exclusion principle. 

It will often be more convenient to label a state not by the quantum numbers n k , 
but rather by h k = n k - (2n)- 1 5^(0) + (2tx)~ 1 5^(0). With these quantum 
numbers on the rhs of ( |2.25p , one should use the phase shifts 5 { k s) {9) of fiTT^ ) on the 
lhs. The advantage of the n k is that now simply n = ^fc=i^/s- We will denote a level 
s(ni, . . . , tin) equivalently as s(hi, . . . , tin)- 

How does one solve eq. ( |2.25|) to obtain the 9 k = 9 k (r) characterizing a given level? In 
general this is of course only possible numerically. One can proceed as follows: For given 
r, type s, and allowed {n k }, rewrite eqs. ( |2.25| ) as 



and then iterate these equations till they converge to a fixed point. One subtlety is that 
when the phases 8^ (6) are determined from a numerical diagonalization of T k (9), one 
must make sure not to "loose track" of a phase when it intersects other phases, or ranges 
over a region larger than 2tt. For the iteration to converge it is perhaps necessary to modify 
the most naive recipe just described on a case by case basis. 

A very simple example, that in practice is quite important and can be treated rather 
explicitly, is that of 2-kink states in the sector of zero total momentum. With n = n\ G 
Z + |, 9 = 9\. 5(9) = 8( a \9\9, —9), we have the following parametric representation of 
the 2-kink energy in finite volume 




(2.26) 





(2.27) 




(2.28) 
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Solving for 6 = 9{r) in a - expansion we find 




2nn 



) 



2 




7r 2 n 2 45 /3 + (47r 2 n 2 /3)(5 



/// 



45') 



= 2 + 



in 



r 



(2.29) 



where 2nh = 2nn — 5, and 5,5', 5"' are the phase shift and its derivatives evaluated at 
9 = (all even derivatives vanish there). This expansion can of course easily be extended to 
higher orders. Even though this expansion has a finite radius of convergence around - = 0, 
it is usually not particularly helpful in practice when comparing with numerical results for 
the finite-volume spectrum obtained using a "small-volume method" like the TCSA, since 
the latter typically becomes inaccurate before the above expansion is useful. Instead, one 
should employ ( |2.28| ) directly, which is correct for all r up to terms exponentially small 
in r; then there is usually a sizable "window of overlap" with the numerical small- volume 
methods, as we will see in the next sections. 

Note that for arbitrary levels the first nontrivial term in the large-r expansion is 
0{r~ 2 ), because of the exclusion principle. Two different levels with the same {n/c} will 
generically be split at 0(r~ 3 ), which always dominates off-shell and tunneling contribu- 
tions. 

Finally, we remark that even though energywise the infinite-volume limit is rather 
trivial — any given level with iV particles approaches Nm — the corresponding eigenstates, 
being eigenvectors of T(0|0), remain nontrivial and distinct, generically. Note that T(0|0) 
is a djsr x permutation matrix; its eigenvalues are therefore some collection of d^-th 
roots of unity. 

3. ^13-Perturbed Unitary Minimal CFTs 

As for any other perturbed CFTs, we have to specify the partition function of the 
unperturbed CFT and the direction of the perturbation. This is crucial for determining 
the nature of the theory, in particular its global symmetry (which may be spontaneously 
broken in infinite volume) and whether the theory is massive or massless. To simplify 
the discussion below we will use the following notation: Let X m ± <p P;q denote the (f)p,q- 
perturbation, in the positive or negative direction respectively, of the unitary minimal 
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CFT [^9| of central charge c m = 1 — m ( T ^ +1 ) and modular invariant partition function 
(MIPF) in the X = A, D, E series [13.0 

Let us first briefly mention the different types of arguments which are helpful in discov- 
ering the kink structure of perturbed CFTs. They give us a better qualitative understand- 
ing of the QFT, and sometimes even provide crucial hints for constructing its S'-matrix. 
(i) Landau- Ginzburg (LG) formulation: It has proven useful in many cases to describe 
CFTs by LG potentials involving a small number of "fundamental" fields |JT| []33| [33 



The perturbed CFT is then obtained by adding to the unperturbed potential a term 
corresponding to the LG realization (in the CFT) of the perturbing field. In the case 
of supersymmetry-preserving integrable perturbations of N=2 super-CFTs (for which the 
potential is believed not to renormalize) , the resulting equations of motion can be treated 
like classical field equations and analyzed for kink solutions. For several models H [[20| the 



results have been used directly to construct the exact quantum spectrum and ^-matrix. 
For theories without N=2 SUSY the LG formulation is of more limited use, since the 
structure of composite operators and the renormalization properties of the theories are 
more complicated. Nevertheless, in some cases it does provide intuition about the kink 
structure [Q. 

(ii) Quantum-group motivated restrictions: In refs. PI ||35|| P| certain perturbed CFTs are 
identified as restricted lagrangian QFTs having SL q (2) symmetry. The restrictions are 
made at specific couplings in the lagrangian which correspond to q being a root of unity. 
The kink structure of the restricted theory is then dictated by the tensor product rules 
of the (finite number of) non-singular finite-dimensional highest- weight irreps of the al- 
gebra U q [su(2)] (see e.g. [|36|]). Namely, the degenerate vacua are labelled by the spins 
j of (possibly a subset of) these representations, and a multiplet of kinks forms a repre- 
sentation of some spin s; kinks in the multiplet interpolate between the vacuum ji and 
vacua j2 appearing in the decomposition of j\ (g> s into U q [su(2)] irreps. The quantum 
group machinary enables one to go far beyond the determination of the kink structure. Its 
disadvantage is that starting from a given perturbed CFT, it is not always straightforward 



10 Since certain <f> Piq are doubled in the D m unperturbed CFTs, one has to specify more precisely 
the perturbation in these cases. The only case this subtlety arises below is that of D5 ± ^1,3, by 
which we will mean the perturbation by the Z2-even combination of the two ^1,3 fields. Also, the 
sign choice for the direction of a perturbation is just a convention for the sign of the perturbing 
field. We use the conventions that are by now standard in the literature. 
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and easy to construct the restricted QFT which it is to be identified with (the intuition for 
choosing the unrestricted theory, coming from the Feigin-Fuchs construction of the CFT, 
is not always clear). 

(iii) Statistical Mechanics: Since perturbed CFTs describe two-dimensional lattice models 
in the scaling region around a critical point, knowledge of the phase structure of the 
off-critical lattice models can give useful hints about the spectrum of the corresponding 
perturbed CFT. 



Returning to the family of 0i ^-perturbed minimal CFTs, we will now briefly review 
some of the known facts about the theories X m ±01,3. As ± 01,3 (= ^-3 ± 02, 1) are 
the two phases of the so-called Ising field theory (IFT) f37j, which can be obtained by 
taking the scaling limit of the Ising model at zero magnetic field from above or below the 
critical temperature, respectively. (From the viewpoint of the underlying lattice model, 
the two theories As ± 0i 5 3 are therefore related by duality.) As + 0i ; 3 is an interacting 
Z2-symmetric theory of a single massive boson with factorizable S- matrix S(9) = — 1. In 
^3 — 0i,3 the Z2 symmetry is spontaneously broken, and the theory is described by the 
simplest (factorizable) ^-matrix theory of a pair of kinks which interpolate between two 
degenerate vacua. Though very simple, it is instructive to discuss this theory in some 
detail, which we will do in subsect. 3.1. 

The theories A m + 0i,3, m > 4, are believed fl38| [[32[ ||3"3|] ||39|1 [pTj to be massless and 



flow to the CFTs A m _i in the infrared (IR). It has recently been argued |[4Q|| [fil] that 
similarly X m + (pi y s flows to X m -i for X = D (m > 6) and X = E (m = 12, 18, 30), while 
D 5 + li3 flows to the CFT A 4 |g88|fl[j. On the other hand, X m - 1)3 is believed to 
be a massive theory for any m > 4 (and any X for which this perturbed CFT is defined). 
Factorizable ^-matrix theories of a multiplet of 2(m — 2) kinks have been proposed for 
A m — 0i,3 independently in 0] (for m = 4, 6) and [|J (for m > 4).0 

In the scattering theories are constructed by a truncation of the Hilbert space of the 
sine-Gordon theory at certain couplings, based on the quantum group (SL q (2)) symmetry 



11 S-matrices for the theories X m — 0i,3 with X ^ A have not been explicitly discussed in the 
literature, except for — <j>i^ for which an S-matrix related through "orbifolding" (cf. J43[ j for a 
statistical mechanics analogy) to that of Aq — (f>\^ was constructed ]IJ. They can presumably [pp| 
be constructed (for m even) by "unitarizing and crossing-symmetrizing" the Boltzmann weights 
of the critical RSOS lattice models of [44]. 
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of the model. The restriction leading to A m — <pi^ is obtained when q = exp(z7r m ? ^ 1 ), 
and so the vacua a are labelled by the spins 0, 1/2, . . . , (ra/2) — 1 of the U q [su{2)] highest 
weight irreps. A kink carries spin 1/2, and so only vacua a and /3 such that \a — (3\ = 1/2 
can be linked by a single kink K a p{9). This situation is graphically described by the 
Dynkin diagram of the simple Lie algebra sl{m) where the dots stand for the m — 1 vacua 
and the links for the pair of kinks interpolating between them. The number of types of 
iV-kink states on a circle, d^, can easily be expressed in terms of the incidence matrix 
I of this diagram: Defining I a p = 1 (0) if the vacua a and (3 are (not) linked, we have 
<ijv = tr(I N ) = Yl^ N i where A are the eigenvalues of /. For the theory A m — ^1,3, 
therefore, d^ = YlT=i [2 cos(7rA;/m)] JV (note that d^ = if N is odd, as it should). 
The S-matrix conjectured for A m — <pi^ is given by M 



U(9) / [2 7 + l] [2(5 + 1] 



e 



a ? k ! 7Ti V [2a + 1] [2/3+1 

X < ' [to + ^tl] ) ' Sinh (^) *<* + Sinh ' °* 



(3.1) 



where 



um = r(I) r(i + J?-) rfi + fl *>W f^-" 

V ; U/ V Tim) V Tim / AJ- # fc (0) i?fc(7rz) 

p/ 2fe 1 zfl \ Y(\ 4- — -U — ) 

v m 7rm / \ m 7rm / 



' 2fc 1 ifl \ 1 2fc 

2fc+l 1 p/i 1 2fc-l 1 fg > 

m -Km' \ m 7rm^ 



(3.2) 



and we use the g-number notation 



a = a 



9 ~ „ _ „-i 



l) a_1 . 7^ (3.3) 



q- 1 v 7 sin(^) 



(which is unambiguous for integer a). The S-matrix (|3.1|) is invariant under the Z2- 
symmetry operation that takes any vacuum a to (m/2) — 1 — a. Note the "crossing 
factors" in the S'-matrix ( j3+l) (cf. the discussion in subsect. 2.3) which ensure strict crossing 
symmetry (I 



Before embarking on a discussion of specific models in the series A m — ^1,3, let us 
point out one feature that is common to the multi-particle spectrum of all these theories. 
Namely, as we now show, the eigenvalues of the iV-particle transfer matrix always come in 
pairs of opposite sign: Note that any iV-kink state, written as ol\ 012 ■ ■ ■ ccn as a shorthand 
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for \K aia2 (9i) . . . K aNOCl (6* at)), has integer (half-integer) either in all even or in all odd 
positions. In other words, there are two classes of vacuum configurations. Let us order 
our basis states such that the first half have integer cti in the even positions i = 2, 4, . . ., 
say, and the second half have integer cti in the odd positions. Due to the restrictions on 
neighboring vacua it is clear that in this basis the transfer matrix is of off-diagonal block 
form ( ^ ^ ) This implies that for any eigenvector of the transfer matrix, written 

in block form as ( Vl | , there is another one f Vl | of opposite eigenvalue, as claimed. 

\v 2 J \~ v 2 J 

When m is odd the two classes of vacuum configurations are actually exchanged by the 

Z2 symmetry of the theory, and so the basis of states can be chosen so that V\ = v 2 are 

eigenvectors of Ai = A 2 ; each pair of eigenstates consists therefore of a Z2-even and a 

Z2-odd state. This is not the case when m is even, where the two eigenstates in each pair 

are of the same parity. 

3.1. Ising Field Theory in the Phase of Spontaneously Broken Symmetry 

This is the case m = 3 of the above family of theories. There are only two vacua, la- 
belled by a = 0, 1/2, which are exchanged by the (spontaneously broken in infinite volume) 
Z2 symmetry. These vacua correspond to the all-spins-up and all-spins-down degenerate 
ground states of the underlying Ising lattice model. There are only two allowed types 
of 2-kink states, 1-^0,1/2(^1)^1/2,0(^2)) and 1-^1/2,0(^1)^0,1/2(^2)), and the two allowed 
scattering processes are given by the amplitudes 5 ,1 /2 ' 1 /2 (6') = Sy£y 2 (0) = -l. 

In general, there are two types of states of any even number iV of kinks on the 
circle, | if 0,1/2-^1/2,0 • • • #1/2,0) and |-#i/2,o#o,i/2 • • • #0,1/2), and none when N is odd. So 
oIn = 2 (0) for iV even (odd). The above two types for iV even transform into each 
other under the Z2 symmetry. Their symmetric and antisymmetric superpositions are the 
eigenvectors of the transfer matrix T = ( J , with eigenvalues 1 and — 1 respectively. 



We therefore have 5^ = and 5^ = ir. The BA equations ( 2.25 ) immediately lead to 



the quantization msinh#fc = 2nhk/R of the momenta in any allowed multi-kink state, 
where the n& are half-odd-integers for the Z2-even sector (type s = 1) and integers for the 
Z2-odd sector (s = 2). The corresponding energies and momenta are therefore simply 



E{ R) = ±,L + ( 2 -f-) 2 , P=ff>, (3.4) 

; — 1 V V / 7 — 1 



k=l 1 x ' k=l 
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where iV is restricted to be even and the are non-coinciding integers or half-odd-integers, 
depending on the sector as above {P is therefore always an integral multiple of 2tc/R). In 
addition to these "peculiar" restrictions, the difference between the energy levels (|3.4|) and 
those of a free particle is that the energy in ( |3.4| ) is not exact for the Ising field theory 
(IFT) when off-shell effects are taken into account. 

The Ising field theory in its high- and low-temperature phases (denoted henceforth 
as iFTM) is to our knowledge the only nontrivial massive QFTs for which the exact 
energy levels are known. They can be obtained from the partition function on a cylinder 
of circumference R and finite length L, with periodic boundary conditions on the spin field 
in both directions, which can be written as |45|] ||15|| 



#)(j;r) = ^ eo(r) { JJ (l + ? e " (r) ) + Y[ (1-Q e " (r) ) 

neZ+l/2 neZ+1/2 

+ (f l{r) [Y[( 1 + 9 e " (r) ) q en(r) )] } 



n£Z n£Z 

Here q = e~ 2irL / R , r = mR (m being the mass of the particle/kink in the +/— case), 



e «( r ) = VW+^' ei(r) = ei(r) - e (r), and 

/oo 
dd cosh6> ln(l±e- rcoshe ) (3.6) 
-oo 

are the scaling functions corresponding to the ground state and the first excitation energies 
Eo,i(R) = (2vr/i?)eo,i(r) in IFT^ - ^ (both decaying exponentially to zero as r — > oo). 

Expanding the expression inside the braces in (|3.5| ) as ^2 Q e ^ r \ one can read off all the 
scaled energy gaps 

^(+\ , . t+\, . , . /x fO Z o -even sector ,„ 

l±) r = e (±) (r) - e (r) = > eJr) + , , * , (3.7) 

w K ' uw ^ nw I ei(r) Z 2 -odd sector . v ' 

Here for the Z2-even sector, which is the same in both phases of the theory, is any 

set of an even number (l-AA^I = 0,2,4, . . .) of distinct half-odd-integers. For the Z 2 -odd 
sector the are sets of distinct integers with |A/" ( - + * ) | = 1, 3, 5, . . . but \Af(~^ | = 0, 2, 4, . . .. 
Comparing these expressions with ( |3.4|) , we see that the latter gives the exact energy gaps 
in the Z2-even sector of IFT^ - - 1 . In the odd sector, on the other the hand, the results 



are not exact. But the exponentially small deviation (2tc / R)ei(mR) = y e mR [1 + 
0((mR)- 2 )] is the same for all levels. This feature is apparently special to the IFT. 
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It is interesting to find the correspondence between the multi-particle states in IFT^, 
specified by the BA quantum numbers hk G and the conformal states of the Ising 



CFT. To do so, recall |46j that the UV limit of the exact scaled energy gaps, e(0), are 



scaling dimensions in the UV CFT. We therefore compare ( |3.7| ) for large r, where it agrees 



with the BA results (3^), with its r — > limit. The only "nontrivial" fact we need is 
ei(0) = 1/8, the rest follows from e n (0) = \n\. Concentrating on the kink phase IFT^~\ 
one concludes that the UV scaling dimensions e^ _ ^(0) of Z2-even states are integers, thus 
corresponding to fields in the conformal families of [0^1 = 1} and [<^i,3 = e] (the dimensions 
of the ancestor primary fields are d± = and d s = 1). More specifically, Z2-even iV-kink 
states in the momentum sector P "come from" [1] ([e]) if f + £P is even (odd). This 
decomposition of the Z2-even sector into two subsectors reflects the fact that [1] ([e]) is 
even (odd) with respect to the Kramers- Wannier duality. On the other hand, the whole 
Z 2 -odd sector originates from [<^i )2 = a] in the UV, as d a = 1/8. 

The characterization of multi-particle states in terms of their quantum numbers hk 
on the one hand and their UV limits on the other hand, leads to the following identities 
for the level degeneracy in c = \ unitary highest weight irreps of the Virasoro algebra: 

Xn^ = (^partitions of n into distinct positive half odd integers) 
Xn^ 2 ^ — (#partitions of n + \ into distinct positive half odd integers) (3-8) 
^(1/16) _ (^partitions of n into distinct positive integers) . 

Here n is a non-negative integer (we set = X^ 16 ^ = -Q) an< ^ X^Hq) = 

q A- 1/48 Xn^\ n is the character of the Virasoro irrep of highest weight A . 



What allows one to obtain the exact partition function ( |3.5[ ), and consequently (|3.8p , 
is the fact that the IFT can be constructed from the theory of a free massive Majorana 
fermion by appropriately summing over the four different boundary conditions for the 
fermion field, cf. eq. (|3.5| ). The case of other integrable perturbed CFTs is much harder; 
the exact (in any analytic form) full finite-volume spectrum is certainly out of reach by 
present methods. However, obtaining the correspondence between multi-particle states 
and UV conformal fields, and ultimately the analog of ( |3.8| ), does not really require the 
knowledge of the exact spectrum at all volumes. The (numerical) TCSA studies together 
with the large- volume BA analysis described in sect. 2 allow us to find this correspondence 
for the low-lying levels, as we will see in the examples studied in the rest of this paper. 
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Our results will extend earlier empirical observations [22] for certain perturbed CFTs 
described by diagonal scattering theories of "ordinary" particles: It was found that in 
certain such cases (though not all) P=0 states of two lightest particles all seem to come 
from a single conformal family in the UV CFT. Furthermore, in these cases the BA results 
rather unexpectedly become exact in the -R — > limit. In IFT( + ) this phenomenon is 
explicitly demonstrated by eqs. ( |3.4p and (^.7|): Zero momentum 2-particle states are in 
the Z2-even sector, the BA equations give their exact energies at any R, and from the more 
general comments above we identify all these states as originating from [e] in the UV. 

3.2. The Subleading Thermal Perturbation of the Tricritical Ising CFT 

We turn now to the theory A4 — ^1,3. Let us label the vacua by a = — , 0, + instead 
of 0, 1/2, 1, respectively, as in (|3.1|) . We also let a, a' take only the values +, — . There are 
eight possible types of 2-kink processes, described by the amplitudes Sqq (9) and S®° a ,(9) 
which can be written explicitly as 

S™'(9) = -2& R(9) W™'{6) , (0) = -2~& R(9) W aa ,{9) , (3.9) 

where 



W ±± (9) = -^=W±±(ni-9) = cosh^ 
W ±T (9) = -^=W± T (ni-9) = -zsinh^ 
and [P]P^[] (see the first reference in for the factors (k — |)) 

«w-n (*-,) r{k _ BnkH+B 

1 I i f°° dx sin(x6'/7r) 

exp < 



(3.10) 



(3.11) 



v/cosh(#/2) "^{4 Jo x cosh 2 (x/2) 
(For explicit computations the alternative form 

is most convenient.) Eq. (|3.11|) is the "minimal" solution to the constraints R(9) = R(ni — 
9) = [cosh{9/2)R{-9)]- 1 . 
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The number of types of iV-kink states in finite volume with periodic boundary con- 
ditions is if iV is o dd, and d N = 2W 2 ) +1 if JV(> 0) is even, of which half are even and 
half odd with respect to the Z2 symmetry which exchanges the a = +, — vacua. 

Now consider 2-kink states of zero momentum, for which we use the basis 

{\K Q .(e)K.o(-9)), \K 0+ (9)K +0 (-e)), \K.q(0)K q .(-0)), \K +0 (9)K 0+ (-9))} (3.13) 

where 9 > 0. In this basis the 2-kink transfer matrix appearing in the first of the two 
Bethe-Yang equations (|2.11] ) is 



T(9\9,-9) = (T(-9\9,-9)y 
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(29) 



The Z2-odd eigenvectors are (1, — 1, ±a/(o — 6)/c, =Fa/(o — &)/c) T , while the Z 2 -even 
ones are (1, 1, ±a/(<2 + b)/c, ±y(a + b)/c) T . The corresponding eigenvalues are given 
by ±y/ (o — b)c and ±a/ (a + b)c, respectively, which explicitly read ±i?(26')-\/cosh6' and 
±z J R(2^) v /cosh^ fi/ 2 {0) where f a (9) = sinh(^f^)/ sinh(^f^). As shown in general 
earlier, the phases come in two pairs, with a difference of 7r between the phases in a pair. 
Explicitly, one of the phases in the Z 2 -odd sector is 



8^ (9) = 5 R (29) , 

while for the Z2-even sector 

5^(9) = 5 R (29) - iarctan(sinh#) 



(3.15) 



(3.16) 



where 5r(9), the phase of R(9) with the branch choice 5r(0) = 0, can be easily read off 
from ( |3.12| ). These phases are shown in fig. 5. Note that 5^ A \±oo) = —8^ B \±oo) = ±-|. 
The parametric form of the zero momentum 2-kink energies is then 



Inn — r)( s U9) 

(r, E) = ( 2mcosh6' ) , se{A,B}, n 6 . 

sinh^ 



(3.17) 



Remember that allowing (positive) integer quantum numbers n as well as half-odd-integers 
in the above is due to the fact that the eigenvalues of (|3.14| ) come in pairs of opposite sign. 
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We now briefly discuss 4-particle levels. There are eight (^-dependent) 4-kink states, 
so we have four pairs of phases, differing by 7r within a pair. Allowing the rik labelling a 
state to be either all integers or half-odd-integers, we only have to consider four types of 
phases, which we will denote by s = A, B, C, D. Types A and C correspond to Z2-even 
states whereas B and D to Z2-odd ones. We will not present analytic expressions for these 
phases, since they are rather messy. Instead, in fig. 6 we show a typical "cross section" of 
the diagonal phases. Note that if we use the quantum numbers hk (cf. subsect. 2.4), then 
all hk G Z or Z + | for type A, £?, C, whereas for type D all hk £ Z + \ or Z — |. In 
particular, in the zero momentum sector all levels of type D will be doubly degenerate. The 
explicit discussion of the lowest 4-particle levels will be presented in the next subsection 
where we compare them with numerical results. 

3.3. Comparison with the TCSA 

The truncated conformal space approach || (TCSA) is a non-perturbative method to 
calculate the low-lying energy levels of a perturbed CFT in finite volume. The idea is to 
truncate the Hilbert space of the CFT to a finite-dimensional subspace on which the hamil- 
tonian of the perturbed theory can be diagonalized numerically. This is to be contrasted 
with conformal perturbation theory (CPT), where the theory is treated perturbatively in 
the coupling constant A, but the Hilbert space is not truncated. The hamiltonian, obtained 
from the perturbed CFT action Acft + A J $ (where the integration is over a cylinder of 
circumference R), simply acts on the truncated Hilbert space as an ordinary matrix 

(j\(H + \V)(R)\i) = ^((A.+A,-^) ^ + (2^ XIV C j9i S Ai _ AiAj _ A .) . (3.18) 

Here are (orthonormalized) conformal states, created by fields of left (right) dimension 
Aj (Aj), y = 2 — (A$ + A$) is the RG eigenvalue of the perturbing field $ which is assumed 
to be relevant and spinless, i.e. y > and A$ = A$, and Cj$>i are the OPE coefficients [^9[ 
of the corresponding fields on the plane. 

As the truncation is removed the eigenvalues Ei(R) of the above matrix should con- 
verge to the exact finite- volume energy levels. There is one subtlety [17] though. If 



y < 1 the (bare) hamiltonian Hq + XV will suffer from UV divergences, so that all of its 
eigenvalues actually diverge as the truncation is removed. If one wants to avoid the ex- 
plicit (and hard) renormalization of these divergences, one can still calculate energy gaps 
Ei(R) = Ei(R) — Eq(R) using the TCSA, since the divergences cancel in the differences 
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of eigenvalues. In view of this fact it is not surprising that even if y > 1 energy gaps seem 
to converge faster than the energy levels as the truncation is removed. We will therefore 
always calculate gaps. Note that y = g for — ^1,3, so this is a case where the energy 
levels themselves will not converge as the truncation is removed. 

The truncation of the Hilbert space is conveniently performed by e.g. ignoring all states 
above a certain scaling dimension or level. Truncation effects become important when 
considering high levels, or if the interaction term gets large compared to the unperturbed 
term in ( |3.18| ). Since the interaction term is proportional to XR y , the TCSA will deteriorate 
with increasing R. 

The perturbed theory is necessarily non-scale-invariant, |A| 1//y being related to the 
mass scale m of the theory. In the cases we will consider, m will be chosen to be the mass 
of the kink-multiplet. The coefficient k relating A to to, |A| = Km y , has been determined 
for some perturbed CFTs to high precision by comparison with the thermodynamic Bethe 
Ansatz (TBA), in particular « = |A|m" 4 / 5 = 0.148695516112(3) for A 4 - </>i, 3 [0. In 
other cases the mass gap has to be determined by estimating the R —>■ oo limit of levels 
corresponding to 1- or 2-particle states. For more details about the TCSA and some case 
studies we refer the reader to P||48||pg|]|]ll||[|r7| (in particular, TCSA results for A4 — 
were first obtained in pBf). 

In fig. 7 we show the first 35 scaled energy gaps ei(r) = -^Ei(R), r = mR, for A4 — (pi^ 
as obtained from the TCSA in the zero momentum sector, together with our BA results 
for 2- and 4-kink levels. Note that at r = the TCSA is trivially exact, the ii(0) being the 
scaling dimensions di = Aj + Aj of fields in the UV CFT, and that the UV fields relevant 
for the 0- momentum sector are all spinless, since up to a factor of 27r/R the momentum 
of a state on the cylinder becomes the Lorentz spin of the corresponding field in the UV 
limit. The analytical BA results were obtained as described in subsect. 2.4. 

The first two gaps correspond to 0-particle states. In infinite volume there are three 
degenerate vacua, but in finite volume they split exponentially, due to tunneling between 



the three wells of the potential. In [|17| we presented strong evidence for an exact expression 
for the first gap ei(r) in terms of the solution of a TBA-like integral equation (cf. also ||16|| ). 
As this is, in particular, an example where the precise form of the large-r behaviour of a 
0-particle level is believed to be known, let us mention that this integral equation predicts 

e x ( r ) = — -L— ifi(r) - (V^ ~ h - - + — + Of-)) . (3.19) 
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Unfortunately, as one sees in the figure, the large-r behaviour is swamped by truncation 
errors. 

The next few levels are the first of an infinite series of 2-kink states, which soon begin 
to overlap with 4-kink states (states of 6 or more kinks appear only at larger energies not 
shown in the figure). The first of each of the "doublets" of 2-kink states correspond to 
type A, the second to £?, with the same quantum number n; BA results are given for the 
n = 2)1>§>2, | levels. The agreement between the TCSA and the BA results is quite 
impressive for the low-lying levels; for small r, large r, and higher levels, deviations should 
be expected due to the limitations of the two methods. In particular, although we have 
not indicated any estimated errors for the TCSA results in the figure, we have checked 
that these errors (as estimated from the variation with truncation level) can explain all 
deviations at large r. Surprisingly, for some levels the BA results are actually exact in the 
r — > limit. We will come back to this point below. 

The first 4-kink level corresponds to |, — |, — |), and the analytical results agree 
quite spectacularly with the TCSA. For the next 4-kink level, -B(§, \, —\, — § ), the moder- 
ate but noticable deviations from the TCSA are mainly due to the (in)accuracy of the latter. 
We identify the next two TCSA levels shown as corresponding to BA type (7(§, |, — |, — §), 
and the degenerate pair l)(|, |, — |, — |), £)(|, |, — |, — |), respectively. That the latter 
level is exactly degenerate is also clear from the perturbed CFT point of view, see below. 
Due to the technical difficulties involved in keeping track of the phases of type C and 
D (cf. the crossings in fig. 6) we have not systematically calculated these energy levels, 
but just computed them for a few selected r values to convince ourselves that the above 
identification is correct. It is actually rather clear from fig. 6 that levels of type B and C 
of the same quantum numbers will be almost degenerate, with level C a bit higher, and 
that a level of type D will be quite a bit higher than the former for large r. 

The next three 4-kink levels are all of type A. The first one corresponds to 
A(2, 1, — 1, — 2), the next one to A(|, |, — |, — |), the last one to the degenerate pair 
A(2, 1, 0, —3), A(3, 0, —1, —2). The BA results shown for the first two of these levels are 
somewhat lower than the TCSA results, for which we blame the latter. (As one can see also 
for 2-kink levels, the general tendency of the TCSA is to give too high estimates for the 
energies as one goes to large r and higher levels). Above these levels there is a "combinato- 
rial explosion" of 4-kink states; unfortunately the TCSA is not accurate enough anymore 
to warrant a further comparison with the BA classification (which can be continued ad 
nauseam). 
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We should briefly comment on the errors of the TCSA results. Experience has shown 
that for large r, and in particular when y < 1, TCSA results can exhibit certain spurious 
features, even qualitatively: In fig. 7, for example, the fact that the second TCSA gap 
increases for r > 3 and the crossing of the first two TCSA 2-particle levels at r « 5.6 
are due to truncation errors. Also, since we generically expect the TCSA to become 
increasingly inaccurate for higher levels, it seems somewhat surprising how good the first 
4-kink level agrees with the BA results in fig. 7; for large r, where we can presumably 
trust the BA results, the agreement is much better than for 2-particle levels of comparable 
energy. We will find a similar feature, even more pronounced, in the example studied in 
sect. 4. It is therefore probably not an accident. Instead, it seems to indicate that the 
volume at which the TCSA begins to deteriorate for a given level is not so much determined 
by the value of the energy at that volume, but rather by the UV scaling dimension of the 
level. This is in fact quite plausible, since the TCSA is non-perturbative in r and the only 
input required for this method are CFT data. 

In table 1 we compare the UV (CFT) and IR (BA) classification of the first 21 levels in 
A4 — <Pi,3- For the fields creating the conformal states in the UV limit we use the following 
notation, for any perturbed CFT. The operators , where n G N and % = 1, . . . , Xn 
(=the degeneracy of level n of the Virasoro irrep of highest weight A and central charge 
c), are certain linear combinations of strings of Virasoro generators L_ ni . ..L_ n with 
ri\ > . . . > n p > 1 and YTj=\ n j = n - The coefficients in these linear combinations depend 
on n, A, c, and the perturbing field They are chosen such that the states \i) = c£} | A) 
form a basis of the irrep of highest weight state |A) at level n in which $ is diagonal, 
namely (7^$^ =0 if i\ ^ %i- (The latter requirement is necessary to ensure that the 
UV states in the table are the A — > limit of energy-eigenstates in the perturbed theory.) 
The operators are defined analogously. 

There are two noteworthy features of the results in table 1. The first is the correspon- 
dence between the mechanisms by which exact degeneracies arise from the conformal and 
BA points of view. From the UV point of view, the degeneracy between states C^^C^^cj) 
and C^^C^^cj) with i\ ^ %2 (where 4> is a spinless primary field of dimensions A = A) 
will be preserved to all orders of conformal perturbation theory (CPT) due to parity in- 
variance, and is therefore exact for all r (this is of course also clear from the TCSA point 
of view). From the BA point of view, on the other hand, levels with quantum numbers 
(n\, -n-2, • • • , tin) 7^ ( — fiNi • • • 7 —Ti2, —fii) will be degenerate, and again because of parity 
invariance this degeneracy will be exact for all r, even when the exponential corrections are 
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taken into account. We can form parity eigenstates by (anti-) symmetrizing over the two 
degenerate BA states, whose UV limits are then identified with superpositions of definite 
parity of the corresponding conformal states (and this is the way in which the UV^TR 
identification of degenerate levels should be made in table 1). Note that from both view- 
points all degeneracies of this kind are two- fold. Higher degeneracies are possible if a 
theory has additional symmetries. 

The second feature is the pattern of differences between the BA ej(r) and the exact 
ones. We will restrict our remarks mainly to 2-kink levels, although qualitatively the same 
picture might be true for states with more kinks. Because the tricritical Ising CFT A4 has 
a finite number (six) of primary fields, the pattern of UV scaling dimensions di repeats 
itself mod 2 (though with varying degeneracies) in any sector of fixed spin. Similarly, for 
2-kink levels in a sector of fixed total momentum the BA e^(0) repeat themselves mod 1 
(for generic scattering theories without the pairing of eigenvalues it would be mod 2). The 
first four 2-kink levels correspond to the fields a' , e', L_iL_icr, L-\L-\e in the UV CFT, 
and the differences between their scaling dimensions and the corresponding BA ej(0) are 
0, Jj, |, Jjj, respectively, cf. table 1. This pattern then repeats itself, presumably forever, 
with certain descendents in the same families as the above fields and BA 2-kink levels 
with higher quantum numbers. All other spinless conformal fields lead to 4- and higher 
multi-kink levels. 

Recalling our discussion (subsect. 3.1) of the BA and the exact energy levels in A 3 — 
</>i,3, one is naturally led to ask if also here the differences between the exact and the BA 2- 
kink levels are just four "universal" functions of r, in the pattern observed above for r = 0. 
(In particular, one might wonder if the BA 2-kink levels with UV limit ei(0) = | (mod 2) 
are in fact exact.) This is not the case, though, as an analysis of eq. ( |3.17|) reveals: Since 
the large-# expansion of the phase shifts (|3.15|) - (|3.16| ) involves powers of e~ e multiplied 
by 9, the scaled BA levels have a small r expansion in powers of r multiplied by powers 
of lnr. Furthermore, the log terms do not cancel in differences of e,i(r) corresponding to 
BA levels of the same type, differing only in their quantum numbers. In contrast, the 
small-r expansion of the exact ii(r), that can be obtained from CPT, is in powers of r y = 
r 4 / 5 , without any log terms. [In theories with diagonal S- matrices, where the scattering 
amplitudes are products of certain universal building blocks (see e.g. [[b|), the phases 
entering the 2-particle BA equations can be expanded in powers of e~ 29 . Consequently, 
the scaled BA energies have a small-r expansion in powers of r 2 , which again contradicts 
CPT predictions, in general. The Ising field theory is an exception: The coincidence of the 
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BA and the exact results for the Z2-even sector in this case is consistent with the above 
remarks because y = 1 and Kramer s-Wannier duality forces all odd terms in the expansion 
in r y to vanish.] 



4. </>2,i-Perturbed Unitary Minimal CFTs 

The theories X m ± 02, 1 (in the notation of sect. 3) appear to be massive for all m. 
In A m ±02,1 with m even the perturbation is "magnetic", manifestly breaking the Z2 
symmetry of A m , and so the sign of the perturbation is immaterial. For m odd, on the 
other hand, the two theories A m ± 02, 1 are different, related to each other by duality. 
Smirnov proposed factorizable S'-matrices for the theories A m — 02, 1, m > 4 (as well as 
for the closely related A m — 0i,2 theories), based on SL g (2)-restrictions of the imaginary- 
coupling A 2 affine Toda field theory. Depending on m, the spectrum contains kinks 
(belonging to one or more multiplets) and possibly bound states of kinks ("breathers"). 
In the case of m = 4 the spectrum consists of a single multiplet of three kinks. For 
this theory, describing the subleading magnetic perturbation of the tricritical Ising CFT, 
Zamolodchikov has independently proposed || a scattering theory with the same spectrum 
and kink structure but different amplitudes. One of the purposes of our paper is to decide 
which (if any) of the two S'-matrices for this theory is correct.^ 

The CFTs D m (m > 5 odd) and E m (m = 11, 17, 29) can also be perturbed by 02, i- 
Smirnov has shown || that by an appropriate "orbifolding" of the S-matrix he proposed for 
A$ — 02,i, one can reconstruct the Z3-symmetric diagonal S-matrix of ||51|| ; this scattering 



theory of an ordinary particle and its antiparticle was previously proposed |5!| for the 
theory D$ + 02, i- Except for this case, we are not aware of any explicit discussion of 
scattering theories for X m ± 02, 1, X 7^ A, in the literature. In general, we think that 
"dual" relations between various particle scattering theories and kink theories, and between 
different kink theories, deserve further investigation. 



12 An attempt to resolve this issue has recently been made in [ |50| by comparing TCSA data 
for finite-size mass corrections with analytical expressions. However, in this case tunneling effects 



have to be taken into account (cf. sect. 1), and in [50] an ad hoc form of the corresponding 
contributions was assumed without justification. Therefore we do not find the analysis of ref. 
satisfactory. 
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4-1. The 02,1 -Perturbation of the Tricritical Ising CFT and its S -Matrix 

From here on we concentrate on the theory A4 — <p2,i- The first information about 



the kink structure of this theory came from the TCSA analysis of [§|§. The large- volume 
spectrum in the zero momentum sector indicated two degenerate vacua — even though 
there is no Z2 symmetry that relates them — and a single 1-particle state of exactly half 
(within the numerical accuracy) the energy of the 2-particle threshold. It was therefore 
concluded that the spectrum of the S'-matrix theory consists of two kinks, interpolating 
between the two vacua, and only one bound state of the kinks ( "propagating only in one of 
the two vacua" ) which is degenerate with them in mass. Alternatively, it is more instructive 
to think about the spectrum as consisting of a single multiplet of three kinks Kqi, Kiq, and 
Kn, so that the incidence matrix describing the linking between the two vacua labelled by 



and 1 is I 



[This is the natural notation from the point of view of the SL q (2) 



structure of the theory, as later proposed by Smirnov || , where the kinks carry spin 1 and 
the vacua are characterized by spins and 1 (q = e l7F / 5 in our conventions embodied in 
eq. (H)).] 

Based on the above spectrum and observations, Zamolodchikov proposed |J to con- 
struct the kink-type S'-matrix of the theory using the Boltzmann weights of the critical 
hard-hexagon lattice model |28|. The motivation comes from the fact that the "spin" vari- 
ables of this IRF model obey exactly the desired restriction imposed on the kinks (namely 
the "spins" can be chosen to take the values and 1, and neighboring "spins" are not 
allowed to be both 0). The nonvanishing Boltzmann weights corresponding to a face of 
the square lattice with the "spins" ct, 7, /3, 5 at its vertices (ordered anti-clockwise starting 
at the upper left vertex) are 
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where u is the spectral parameter. Explicitly, this is 
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These Boltzmann weights satisfy the Yang-Baxter equation. 

Zamolodchikov's proposal is to look for a kink ^-matrix of the form 



PiPs 
P a pp 



R{9) W 



a 5 
.7 P. 



(\e) 



(4.3) 



with A, p a , and R{9) suitably chosen to satisfy unitarity 



crossing symmetry ( |2.3|) 



and lead to a direct- channel simple pole corresponding to one of the kinks at 9 



2-wi 
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all scattering amplitudes, except Sqq(9) (since there is no kink K 00 ). 

Actually, a stronger version of this "pole constraint" must hold, namely the bootstrap 
equations ( |2.5|) should be satisfied. It is convenient to postpone the discussion of the 
bootstrap, since already the weaker version of the pole constraint together with crossing 
symmetry can be satisfied only if p a = s a (ignoring an irrelevant ct-independent overall 
factor), and (i) R(in - 9) = +R(6) and A G z(§ + 6Z) or (ii) R(in - 9) = -R(6) and 
A G i(— | + 6Z) .HH Unitarity implies the constraint 
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Zamolodchikov chooses A = — In this case the constraints on R(9) have the "min- 



imal solution" 
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taken in §. Here f a (6) = sinh(^±f^)/ sinh(^f^). 

But A = — 6z/5 is not the only choice. Choosing (with hindsight) A = 9z/5, the 
"simplest" solution for R{9) is 
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A similar ambiguity in the analog of Zamolodchikov's A shows up in Smirnov's quantum- 
group approach. Since he is working with the A\ ' affine Toda field theory for arbitrary imaginary 
coupling, he can fix it as follows: For generic imaginary coupling the model contains a (lightest) 
breather, a kink-kink bound state which is an "ordinary" particle. Smirnov demands that its 
scattering amplitude with itself be the analytic continuation of that of the single boson in the 
real-coupling ASp Toda theory of ref. |53|| . (It is known that the analogous analytic continuation 
indeed relates the scattering amplitude of the single boson in the sinh-Gordon model to that of 
the first breather in the sine-Gordon model.) 
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However, the resulting ^-matrix ( }4.3| ) exhibits some unwanted simple poles in the physical 
strip, namely at 9 = and 9 = There is a cure, though: One can always multiply 
the above R(9) by arbitrary "CDD factors" |2|] of the form Ela^^)' where F a (9) = 



—f a (9)f a (Tvi — 9), which do not affect the algebraic constraints on R{9). Muliplying by 
F_i/g(9) cancels the unwanted poles, and a further factor of F 2 /q(9) cancels the zeros of 
/_2/5(96*/5) and /3/5(96'/5) in the physical strip without introducing new unwanted poles. 
Choosing the '+' sign in ( [4.6| ) to get the right signs for the residues of the poles in the full 
S'-matrix, cf. (|2.6| ), we see that all physical constraints are satisfied by 



in a "minimal" way, in the sense that (|4.7| ) has the smallest number of poles and zeros 
in the physical strip. And this solution corresponds to Smirnov's proposal!0 Smirnov 
arrived at his proposal from a completely different direction || , and we — being aware of 
his result — just rederived it using considerations employed by Zamolodchikov in order to 
clarify the relation between the two approaches. 

The function Rs(9) looks somewhat complicated, in particular since we had to include 
the factors F_ 1 / g (9) and F 2 /g(9). The presence of these factors becomes less mysterious 
when considering the bootstrap equations ( |2.5|) . One can check that they indeed hold for 
Smirnov's S'-matrix (with or without the crossing factors), using 



fls(9 + f)fls( ,_ f) .__^_. (4 , 



Given the factor F_i/ 9 (9), necessary to eliminate the unwanted poles, it is crucial to also 
include the factor ^2/9(6*) in Rs(9), otherwise the bootstrap equations are not satisfied. 

There is one more important point we have to discuss in connection with Smirnov's 
S'-matrix. In ref. |J there are no crossing factors, i.e. all p a = 1 in ( |4.3[ ). Smirnov claims 
that the resulting violation of crossing symmetry does not make the theory inconsistent. 
However, he has not explained how such a violation could arise in a QFT, which is indeed 
rather difficult to see. It would be nice, to perform a direct check of the presence or absence 
of the crossing factors. Unfortunately, as we saw in sect. 2.3, these factors do not affect 



14 In || Rs(0) was given in an integral representation; it was rewritten in (essentially) the 
form (14.71) in ref. 
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the finite- volume spectrum of a QFT (modulo, possibly, exponentially small terms), so 
comparison of BA with TCSA results will not reveal if they are present or not. 

Returning to Zamolodchikov's S-matrix, we note that for Rz{9) eq. ( |4.8| ) holds with 
a plus sign on the rhs, which leads to a violation of all the bootstrap equations ( |2.5| ) by a 
sign. Since it is just one overall sign, we can restore the bootstrap equations by reversing 
the sign of all S-matrix elements. (It would have been a serious problem if some bootstrap 
equations were satisfied and others violated by a sign, because unitarity, crossing and the 
Yang-Baxter equations fix all relative signs.) The overall sign that Zamolodchikov chose || 
was presumably motivated by the fact that it leads to positive (imaginary) residues for the 
direct-channel poles of S'-matrix elements of the form S^{6). Such residues correspond to 
real couplings g ai p in ( |2.6|) , as expected for a unitary theory. On the other hand, his sign 
is such that S]Z(0) = +1 for all allowed a,P,j, which according to ( j2.10| ) means that all 



kinks are fermions. But then it is hard to see how they can be bound states of each other. 
Changing the overall sign of Zamolodchikov's S'-matrix solves this problem (since then 
all kinks are bosons), in addition to restoring the bootstrap equations. The S-matrix will 
however now violate 1-particle unitarity, i.e. some couplings g ai (3 will be purely imaginary. 



According to previous experience [[54] this suggests that Zamolodchikov's S-matrix — if it 
describes any consistent theory — might be that of a perturbed reon-unitary CFT. 

Still, since the a priori problems with Zamolodchikov's S-matrix are due to signs, a 
somewhat subtle issue in S-matrix theory (cf. ||25|| ||14|| [[I If), and since problems with one 
conjecture do not prove that another (Smirnov's) is correct, it is important to provide a 
more direct check of these proposals. This will be done in the next subsections. 

4-2. Multi-kink states 

Imposing periodic boundary conditions the number of different types of A-kink states 
is d N = tr Q\) N = ( i ^) Ar + ( i ^) JV , or recursively d N = d N - X + d N - 2 with d = 2 
and d± = 1. Note that states with an odd number of kinks are allowed here, contrary to 
the theories A m — of sect. 3. In particular, there is one type of 1-kink states |ATii(#)), 
with 6 quantized in finite volume R according to msinh# = 2-Kn/R with h G Z. Up to 
exponentially small corrections the corresponding energy is just y^m 2 + (2nh/R) 2 . 

For 2-kink states of zero momentum on a circle we use the basis 

{ \K Q1 (6)K 10 (-6)), \K 10 (6)K Q1 (-6)), {Kn^K^-d)) } . (4.9) 
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The notation 



a V>) 


— a{V) — dn^) 






KB) 


= pi /[27Ti) b{0) = 


S°il(0) 


= S$(6) 


c(0) 


= Pi c{9) -- 


= S l l(9) 


= Slh(0) 


d{9) 


= p- ej ^ l) d{9) = 


- S&(0) 




e(0) 


= p [/M e{6) = 


Sff(0) 





(4.10) 



(4.12) 



will be useful below. Here p\ = 2 cos ^ for Zamolodchikov's or the crossing-symmetrized 

Smirnov proposal, whereas p\ = 1 for Smirnov's original proposal. Note that the tilded 

quantities would be the S'-matrix elements if there where no crossing factors. 

The 2-kink transfer matrix now reads 

/ e(20) 6(20) \ 
T(0|0,-0) = - d(29) . (4.11) 

\ 6(20) a(20) / 

Its three eigenvalues 

A(0) = -d{29) 
A ± (0) = \(d-~a±\J{d-~a) 2 + 4gj)(20) , 
correspond to the eigenvectors 

v{9) = (p?'™ p[^\ pT f , 
v±(9) = (1/VxI, y/x±, v^^±)(0) T . 

Here 

X±{9) = ~>t(f) ' = - A ± (0) + a(20) • (414) 

Note that in agreement with our general proof in sect. 2.3 the eigenvalues are independent 
of the crossing factors, whereas the eigenvectors are not. 

Let us denote the phases of A(0), A±(0) by <5 (j4) (0), 5 (B) (0), 5 (c) (0), respectively. 
Fig. 8 shows a plot of the three types of phase shifts for Smirnov's ^-matrix. In terms of 
these phases we have the standard parametric form ( |2.28| ) of the 2-kink energy levels. 

Let us now turn to 3-kink levels. For the four-dimensional space of 3- kink states on a 
circle we use the basis 

{ \K 01 (9 1 )K 11 (9 2 )K 10 (9 3 )), \K 11 (9 1 )K 10 (9 2 )K 01 (9 3 )), 

(4.15) 

|ifio(0i)#oi(02)^ii(0 3 )>, |ifii(0i)^ii(02)#ii(0 3 )> } ■ 



(4.13) 
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It is possible to give explicit analytic expression for the eigenvalue phases of the 3-kink 
transfer matrix for general 61,62,03. However, except for one of the eigenvalues (see below) 
they seem to be rather complicated; we will not write them down, since in the end the BA 
equations have to be solved numerically anyhow. Instead, we show in fig. 9 a representative 
"cross section" of the four phase functions for Smirnov's S'-matrix. 

For a special class of 3-kink states the relevant phases have very simple expressions, as 
we will presently discuss. For reference, let us first write down the transfer matrix relevant 
for general 3-kink states: 



T(6\6 1 ,6 2 ,6 3 ) 



( 


cce 


dbb 


c 


a 


b 


bbd 





c e c 


a 


b 


c 


e c c 


bdb 





b 


c 


a 


\ b a c 


a c b 


c b a 


a 


a 


a 



(4.16) 



/ 



where the i-th factor in each of the triplets of functions is to be evaluated at 6 — 6i. Note 
that for 6 G {#1, 62,63} the structure of this matrix simplifies since 6(0) = 0. 

One would expect that in the zero momentum sector there are solutions of the Bethe- 
Yang equations with 62 = 0, and, in fact, that the lowest 3-kink levels are of this form. 
This is indeed true. 6 2 = means that we have to find eigenvectors of T(0\9, 0, —6) with 
eigenvalue equal to —1, where we have set 6 = 9\ = —63. This eigenvalue turns out to have 
a two-dimensional eigenspace for all 6, which we can write (projectively) as (z, 1, 1, (zc(6) — 
e{9)) /b{9)) in terms of a free parameter z. Since transfer matrices T(6) with different 6 
can be simultaneously diagonalized, there must be (at least) two values of z where this 
vector is also an eigenvector of T(9\9, 0, —9). We find that these values are 

corresponding to the eigenvalues 

a +ot = - m - m^i 

KQ) (4.18) 

A_(0) = -c(0)c(20) 

of T(9\9,0,—9). These eigenvalues correspond to the 9 2 = section of the phase shift 
functions labelled A and B, respectively, in fig. 9.0 



15 And this is perhaps the appropriate moment to remark that for general 6i the eigenvalue of 
type B of T fc (0i, 02,0s) is A^ B) (0i, 2 , 3 ) = -c(0 k i)c(9 kj ), where k,i,j £ {1,2,3} are all distinct. 
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In terms of the corresponding phases 5( s \9), ( |2.19| ), the parametric form for the 3-kink 
energies of type A, B with 6*2 = is 

/ 27m — 5^(9) \ 
(r,E) = ( V ; , m (l + 2coshfl) ) , sG{A,B}, n G N . (4.19) 

All other zero momentum 3-kink levels, for which we do not give explicit expressions, 
come in degenerate pairs. They are either the pairs s(hi, 77.2,773) and s(— 77,3, — 77,2, — hi) 
with s = A, B and the hi distinct nonzero integers that sum up to zero, or C(hi, 77,2,773) 
and D{— 77,3, — n 2 , — hi) with distinct rlj G Z + | whose sum vanishes. 

4-3. Comparison with the TCSA 

For A4 — 02,i the coefficient /? = |A|m -9 / 8 has not yet been determined to high 
precision by comparing CPT with the TBA. Performing the TCSA with A = ^- we estimate 
m = 0.99(1) (cf. [p0[ ) by looking at the variation (with truncation level) of the 1-kink energy 
for moderate to large volume. Hence k = 0.161(2). Compared to the case of A4 — ^3 we 
will see that the greater inaccuracy in k for A4 — <p2,i will be more than compensated by 
the larger value of y, so that the TCSA is accurate up to larger values of r. 

In fig. 10 we show the first 21 scaled energy gaps in the zero momentum sector obtained 
from the TCSA for the </>2,i-perturbed tricritical Ising CFT, together with our BA results 
for 2- and 3-kinks as derived from Smirnov's S'-matrix. What is clear right away is that 
the latter S'-matrix receives very strong support from our results; the S'-matrix of gives 
quite a different picture, definitely not describing A4 — 4>2,\- 

Since there are two degenerate vacua in infinite volume, there is now only one gap 
approaching zero energy exponentially. The next gap corresponds to the kink K\\ at zero 
momentum, which is followed by a series of 2- and 3-kink states, and we can also see the 
first three 4-kink levels (for which we have not performed a BA analysis, mainly out of 
laziness). The UV and IR classification of levels is shown in table 2. 

Because of our detailed discussion of the A4 — 0x,3 case we can be more brief here. 
Now the differences between the 2-kink BA e^(r) and the exact gaps at r = are Jj, 0, Jj, 
a pattern which is repeated mod 2. So it appears that 2-kink states of type A, £?, C 
correspond to fields in the families [a'], [e'}, [a], respectively, of the tricritical Ising CFT. 
One can similarly speculate how the 3-kink states of different types correspond to specific 
conformal families: The 3-kink levels of types A and B with #2 = correspond to the 
families [e] and [e"\ , respectively, whereas the "conjugate" types C and D are associated to 
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both of the families [a] and [e'}. (It is not so clear, though, what the UV-IR correspondence 
is for 3-kink levels of type A and B with 62 7^ 0.) Again, as in subsect. 3.3, one can show 
that the differences between the exact and the BA energy levels can not be some finite 
number of "universal" functions (corresponding to the types of BA levels, independent of 
the {nk}), as they are in the Ising field theory of sect. 3.1. 



5. Discussion 

The multi-particle finite- volume spectrum provides a very characteristic "fingerprint" 
of a QFT. We discussed how to calculate this spectrum for large volume in integrable 
scattering theories of kinks. For simplicity we considered theories with a particle spectrum 
consisting of a single multiplet of kinks, in finite volume with periodic boundary conditions; 
these assumptions can easily be generalized. As comparison with numerical results from 
the truncated conformal space approach has shown in several cases, the analytical large- 
volume predictions are usually even better than one could have hoped for, in that they are 
accurate up to quite small volumes (in fact, for some levels they become exact again at 
zero volume). 

Applying our results to the case of the subleading magnetic perturbation of the tricrit- 
ical Ising CFT, comparison with TCSA data provides very strong support for the S'-matrix 
conjectured by Smirnov || ; the S'-matrix conjectured by Zamolodchikov || definitely does 
not describe this theory (as suggested in subsect. 4.1, it might describe a perturbation of a 
non-unitary CFT). There remains, however, one subtle question that cannot be answered 
by this comparison. Namely, if the original S'-matrix of Smirnov should be multiplied by 
"crossing factors" in order to restore crossing symmetry in the "traditional" form ( |2.3|) . 
These crossing factors do not affect the multi-particle spectrum (at least not the dominant 
terms). The same is true for 1-particle levels, and the exact ground state energy which 



can be calculated ]55| using the thermodynamic Bethe Ansatz is easily seen not to depend 
on the crossing factors for any volume. Nevertheless, the S'-matrices are different and can 
in fact be distinguished "experimentally" , by measuring, for example, time delays (which 
are proportional to phase shift derivatives) in scattering experiments. Since it is hard to 
see how a non-crossing-symmetric S'-matrix could correspond to a consistent QFT, one 
might favor the crossing-symmetrized S'-matrix as the correct one. Still, it would be nice 
to perform some explicit check. 
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Shifting to more general issues, an important open problem in the study of (massive) 
integrable perturbations of CFTs is the relation between the UV CFT classification of 
states and the IR Bethe Ansatz classification. For the examples we have studied, our 
results show explicitly how this "UV to IR map" works for the lowest 20 levels or so in 
the zero momentum sector. In addition to extending this map to arbitrary levels and 
figuring out the combinatorics involved (will this lead to new ways of writing Virasoro 
characters?), one wonders if it is possible to understand more conceptually why a given 
multi-particle state corresponds to a certain conformal field in the UV. Finally, it remains 
to be seen if the simple universal pattern that holds for the differences between the exact 
and the Bethe Ansatz finite-volume levels in the Ising field theory generalizes to a more 
sophisticated pattern in other integrable QFTs. 

Throughout the paper we considered only integrable theories. At least for 2-particle 
states this restriction is not really necessary. By increasing the volume the energy of any 
level in a generic non-integrable theory will eventually drop below the threshold for parti- 
cle creation, with only 2-particle states scattering elastically among themselves remaining. 
The results of sect. 2 then apply without any modification. Of course, for a non-integrable 
theory there is not much hope of ever knowing the exact S'-matrix. Therefore one might 
want to use our results "in reverse" , namely extract the scattering amplitudes by compari- 
son with numerical data for the energy levels |H| f2^] |2^] (and this is just as interesting for 
integrable theories whose ^-matrix is not known). Looking at 2-particle levels in the zero 
momentum sector, one can directly only obtain the diagonal phases 6^ (9) = 8^ a \d\9, —9). 
Unfortunately, if one is just given the energy levels, one does not know how to relate the 
diagonal to the asymptotic 2-particle basis of states. Except in rather special cases, where 
the theory is a priori known to have a large (global) symmetry so that the relation between 
these bases is very simple, it will be quite difficult if not impossible to extract the physical 
scattering amplitudes from the diagonal phases 5^ (9). 

We mentioned in the introduction that the sine-Gordon (SG) model contains solitons 
which are not kinks, in that there are no nontrivial restrictions on the multi-particle Hilbert 
space. We would like to conclude with some remarks on certain variants of the SG model 
that do contain kinks. Note that to define the (classical) SG theory we do not only have 
to specify its lagrangian, ^d^tpd^tp — (mo//3) 2 (l — cos (3(f), but also the "target space" 
on which the field (p lives. For the standard SG model (p lives on a circle of radius r' = 4, 
i.e. one identifies field configurations differing by a period of the potential. But what if 
we let (p live on a circle of radius r' = ^, for some generic Let us denote the 
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This possibility was also noticed by Swieca [56|. 
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corresponding theory by SG((3, k). All properties of the ordinary SG model SG((3, 1) that 
rely only on local properties of the field (p, like classical integrability, will also hold for 
SG((3,k). We expect similar statements to be true also at the quantum level (provided 
that P 2 < 8tt). In a semi-classical framework it is clear that the quantum theory SG(f3, k) 
will have a /c-fold degenerate vacuum. Labelling the vacua by a = 0, i, . . . , the theory 
will have kinks K aa > with \a — a!\ = \ (or ^^). As in the ordinary SG model there may 
also be bound states of kinks. 

The S'-matrices of SG(@, k) with different k are also closely related: Using the ob- 
vious Zfc symmetry, as well as time-reversal and parity invariance, we see that every 
non-vanishing 2-kink amplitude is equal to one of the three (a-independent) amplitudes 
S2££(0), Sa + a Q+ (#), Sa + a Q "(#), where a± = a ± §, a ++ = a + 1 (mod *fi). These 
should be independent of the global properties of the field (p and therefore equal to the 
soliton-soliton, anti-soliton-soliton reflection, and anti-soliton-soliton transmission ampli- 
tudes, respectively, of the standard SG model fl24j] . Nevertheless, the global properties of 



if, reflected in the fact that there are restrictions on multi-particle states, do give rise to 
differences between the theories with different k, e.g. in their finite- volume spectra, as is 
clear by just contemplating the analysis of sect. 2 for different k. 

Finally, it is illuminating to consider the UV limit of the theories SG(f3, k). Changing 
to standard CFT conventions, where the massless limit of the (euclidean) SG field y/ntp is 
denoted by X, we see that the UV limit of SG(j3, k) is the c = 1 gaussian CFT at radius 
r = (in the notation of |p7|| )Jl3 Gaussian models at different r differ in their field 

content, and considering SG(f3, k) as a perturbed CFT, namely the perturbation by the 
dimension \{^:) 2 = j- spinless vertex operator cos(kX/r), it is again clear that theories 
with different k will have different finite-volume spectra (although the ground state energy 
is the same) . We plan to provide a more detailed analysis of various aspects of the theories 
SG(P,k) elsewhere. 
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17 Note that all these CFTs can be obtained from the one with r = \J~f3/i{ by "orbifolding" 
with respect to the Zfe subgroup of the U(l) ("winding") symmetry of the gaussian model. By 

analogy, SG({3, k) can be thought of as a Zfe "massive orbifold" of SG{(5, 1). 
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Table 1: UV and IR classification of levels in A4 — 0i,3- 
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1(1,0,-1) 


2 


L_ 1 L_ 1 a' 


2+1 


^d' — f ) 


2+! 


e" = 03,i 


3 


5(1,0,-1) 


2+f 


L_iL_i e' 


2+ f 




2+f 


L—2L—2 1 


4 


1st 4-kink state 




r (3/80) ~(3/80) 
/ -2,1 ^2,2 

^(3/80)^(3/80) 
^2,2 ^2,1 


4 + — 

^ ^ 40 

4 + — 

* ^ 40 


/V4 1 5\ 
°V3' 3' 3' 

V3' 3' 3^ 


4 
4 


r (3/80) ~(3/80) 
^2,1 ^2,1 


4 + — 

^ ^ 40 


(7(2,-2) 


4 


r (3/80) ~(3/80) 
^2,2 ^2,2 


4 + — 

^ ^ 40 


2nd 4-kink state 




L_2-L_2 £ 


4+1 


1(2,0,-2) 


4 


L-2L-2 a' 


4+1 


^(f> _ |) 


4+f 


L-xL-x e" 


2+3 


5(2,0,-2) 


4+f 


-(3/5)^(3/5) / 

^2,1 ^2,2 fc 

,.(3/5) =(3/5) / 
^2,2 ^2,1 fc 


4+ f 
4+f 


£>(5 2 _7x 

V3' 3' 3' 

^( 3' _ 3' _ 3) 


4+f 
4+f 


r (3/5) *(3/5) / 
^2,1 ^2,1 e 


4+f 




4+f 


.(3/5)^(3/5) , 
^2,2 ^2,2 fc 


4+f 


3rd 4-kink state 





Table 2: UV and IR classification of levels in A 4 — <fi 2 ,i- 
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Figure Captions 



Fig. 5: Diagonal 2-kink phases S^ s \9) = 5^ s \9\9, —9) (in units of n) for the S-matrix 
(^9| ) - (|3.11|) , conjectured for A 4 — (j>\^. Shown are type A, eq. ( |3.15| ), as the upper line, 
and type £?, eq. ( |3.16|) , as the lower line. 

Fig. 6: The "cross section" i— »■ ^S^(9\9, \9, —\Q, —9) of the diagonal 4-kink phases of 
type s = A, B, C, D (labelled from top to bottom for small 9) for the S-matrix ( |3.9| )-( |3TTl| ). 

Fig. 7: Scaled energy gaps e^(r) in the zero momentum sector as obtained from the TCSA 
for the 0i 5 3-perturbed tricritical Ising CFT and the BA using the S'-matrix (|3.9|) -( |3?TT| ). 
Shown are the first 35 gaps calculated with the TCSA using the 228 states in the CFT 
up to level 5 (solid lines), the first ten 2-kink levels from the BA (dotted lines), as well as 
BA results for the first 4-kink level of type B (short dashed line) and the first three 4-kink 
levels of type A (dashed lines). The next 4-kink level, for which no BA results are shown, 
is doubly degenerate. 

Fig. 8: The three types s = A, B,C (from top to bottom) of diagonal 2-kink phases 
5( s )(9) = 5( s \9\9, -9) (in units of tt) for Smirnov's S'-matrix conjectured to describe the 
<^2,i-perturbed tricritical Ising CFT. 

Fig. 9: The cross section 9 i-> ^ s \9\9, |0,0 3 (0)) wh ere sinh0 3 (0) = -(sinh0 + sinh ±0), 
for the four types s = C,D, A,B (from top to bottom) of diagonal 3-kink phases in A 4 — <p 2 ,i 
(using Smirnov's S-matrix). 

Fig. 10: Scaled energy gaps e^(r) in the zero momentum sector as obtained from the 
TCSA for the 02, l -perturbed tricritical Ising model and the BA based on Smirnov's S'- 
matrix. In addition to the first 21 TCSA gaps calculated with 228 states (solid lines), we 
show the first eight 2-kink levels from the BA (dotted), and the BA results for the first 
four 3-kink levels with 02 = (dot-dashed), which are alternatingly of type A and B, as 
well as the first two degenerate pairs of 3-kink levels of type C and D (dashed). The levels 
for which no BA results are shown are the first three 4-kink levels. 
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